Power-law cosmic expansion in f{R) gravity models 
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We show that within the class of f{R) gravity theories, FLRW power-law perfect fluid solutions 
only exist for i?" gravity. This significantly restricts the set of exact cosmological solutions which 
have similar properties to what is found in standard General Relativity. 
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INTRODUCTION 

Currently, one of the most popular alternatives to 
the Concordance Model is based on modifications of the 
Einstein-Hilbert action. Such models first became pop- 
ular in the 1980's because it was shown that they natu- 
rally admit a phase of accelerated expansion which could 
be associated with an early universe inflationary phase 
The fact that the phenomenology of Dark Energy 
requires the presence of a similar phase (although only 
a late time one) has recently revived interest in these 
models. In particular, the idea that Dark Energy may 
have a geometrical origin, i.e., that there is a connection 
between Dark Energy and a non-standard behavior of 
gravitation on cosmological scales is consequently a very 
active area of research. 

One such modification is based on gravitational actions 
which are non-linear in the Ricci curvature R and/or 
contain terms involving combinations of derivatives of 
R [3, 0, 0, [^. Over the past few years, these theories 
have provided a number of very interesting results on 
both cosmological [1, 0, [1, H, and astrophysical d, 11 1 
scales. An important feature of these theories is that the 
field equations can be recast in a way that the higher or- 
der corrections are written as an energy - momentum ten- 
sor of geometrical origin describing an "effective" source 
term on the right hand side of the standard Einstein 
field equations In this Curvature Quintessence 

scenario, the cosmic acceleration can be shown to result 
from such a new geometrical contribution to the cosmic 
energy density budget, due to higher order corrections of 
the Hilbert-Einstein Lagrangian. 

Of considerable importance to the study of the cos- 



mology of these models is the existence of exact power- 
law solutions corresponding to phases of cosmic evolu- 
tion when the energy density is dominated by a perfect 
fluid. The existence of such solutions is particularly rele- 
vant because in Friedmann-Lemaitre-Robertson- Walker 
(FLRW) backgrounds, they typically represent asymp- 
totic or intermediate states in the full phase-space of the 
dynamical system representing all possible cosmological 
evolutions. 

In this paper we investigate the implications for the 
gravitational action if exact FLRW power-law solutions 
in f{R) gravity are assumed to exist. We discover that 
such solutions only occur for a very special class of f{R) 
theories. This result is complementary to one recently 
found for /(G) gravity models [li| . 



FIELD EQUATIONS FOR HOMOGENEOUS AND 
ISOTROPIC f{R) MODELS 

We consider the following action within the context of 
four-dimensional homogeneous and isotropic spacetimes, 
i.e., the (FLRW) universes with negligible spatial curva- 
ture: 



A = 



^[f{R) + C. 



(1) 



where R is the Ricci scalar, / is general differentiable (at 
least C^) function of the Ricci scalar and Cm corresponds 
to the matter Lagrangian. Units are chosen so that c — 
16ttG= 1. 

It follows that the held equations for homogeneous and 
isotropic spacetimes are the Raychaudhuri equation 
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p + 3P + f-fR + ef'R + 3f"'R^ + 3f"R ] 



(2) 



where Q is the volume expansion, which deflnes the scale factor a{t) along the fluid flow lines via the standard 
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relation Q = 3a/a, and /^"^ abbreviates d"f /{OR)" for 
n = 1..3; the Friedmann equation 



ef'R 



(3) 



the trace equation 

3Rf" = p-3P + f'R - 2/ - 3Qf"R - 3f"'R^ ; (4) 
and the energy conservation equation for standard matter 

p = -e{p + p). (5) 

Combining the Friedmann and Raychaudhuri equations, 
we obtain 



where m > is a fixed real number. We further as- 
sume that the standard matter can be described by a 
barotropic perfect fluid such that P — wp with w G 
[—1, 1]. From the energy conservation equation, we ob- 
tain 



(8) 



and from (El) we see that the Ricci scalar becomes 



R = 6m(2m — l)t ^ = a-mt 



(9) 



Requirements for the existence of power— law 
solutions 



Analogously to , let us now assume there exists an 
exact power-law solution to the field equations, i.e., the 
scale factor behaves as 



i{t) = aot" 



(6) Note that i? > if m > 1/2, and i? < for < m < 1/2, 
so the value of m fixes the sign of the Ricci scalar. 

Using the background solutions above, we can write 
the Friedmann, Raychaudhuri and trace equations in 
terms of functions of time t only, assuming with no loss 
of generality that t > 0. 

Considering values of m 7^ 1/2, we can then solve ^ 
for t and re-write these equations in terms of the Ricci 
scalar R, f{R) and its derivatives with respect to R. The 

(7) Friedmann equation for example becomes 



f'R' 



-f'R 



1 - 2m 



/-f (2m- l)K 



R 

am 



= 0, 



(10) 



r 



Note that for the power-law so- 



where K 

lution ([7]), R/am is positive at all times, and therefore 
equation (fTO]) is real-valued over the range of R. 

Since we want ([71) to be a solution at all times, i.e.. 



J 



R spans over an entire branch of the real axis, we can 
interpret (jlOp as a differential equation for the function 
/ in R space. Solving this equation gives the following 
general solution 



f{R) - ( — 



2 3 



(11) 



r 



where we have abbreviated 

4(2m- l)po 



A — 



2 - m(13 -t- 9w) + 3m2(4 + 7w + 3^2) 



, (12) 



= 1 + 10m + m^ (13) 
and C1.2 are arbitrary constants of integration. We note 



that the above form of / identically satisfies the other 
field equations, if we similarly convert them into differ- 
ential equations in R space. Note that j3m > for cosmo- 
logically viable solutions with m > 0. This means that 
the exponents in the solution are all real valued in the 
case considered here. Also Amw is real- valued and non- 
zero unless m = 1/2, but diverges if m and w satisfy the 
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relationship it; = (3 — 7m ± y/ Pm) /6rn. In general, the 
function /(i?) is real-valued if m and w do not satisfy the 
above relationship, and if i? > (i.e., m > 1/2). Fur- 
thermore, If we want to ensure that for m = 2/[3(l + w)] 
and K = 4/[3(l + w)^] the theory reduces to GR, then 
we have to set Ci = C2 = 0. In that case, the solution 
(fTTj) is real- valued for all R provided R/am > 0. If we 
re-write |m(l + w) = n, then we can see that we recover 
the well-known result that in i?"-gravity, there exists an 
exact Friedman-like power-law solution a oc t2n/(3(i+-i«)) 
The GR-limit can now be identified as the case n = 1. 

Scalar field analogy 

It is interesting to use the solutions found above to 
reconstruct the effective scalar field often invoked to de- 
scribe the dynamics of f{R) grav ity models. Using this 
analogy, it has been argued in [13] that f{R) theories suf- 
fer from a singularity problem, namely that in the past, 
at finite time, the dynamics drives the model towards in- 
finite values of the curvature corresponding to points in 
the scalar field potential atteignable for finite values of 
the scalar field. Moreover, the effective potential of the 
models studied in [l3| are multivalued, which is a very 
unnatural feature. In what follows we will show that 
the models dTTI) that lead to power-law solutions for the 
scale factor do not suffer from such pathological behav- 
iors, but admit a well-defined scalar field representation 
with a single- valued potential and no curvature singular- 
ity. 

The fact that the curvature is well behaved can be di- 
rectly inferred from Eq. ^ since the only divergence 
occurs for t — 0, or equivalently, for a — 0, and this sim- 
ply corresponds to a standard Big-Bang type singularity. 

We adopt the representation in terms of a scalar field 
used in 13], by defining the scalar field (f> and its potential 
V{(p) through the following equations: 



dR ^ 3 [^^^''^ -dR^'jdR^- 

The shape of the potential is illustrated in Figs. [1] and 
[2] for various values of the equation of state and of the 
constants Ci and C2. 




FIG. 1: Upper left plot: 1/(0) for «; = and 3m/2 = 1.1 
in the case Ci = C2 = 0. This corresponds to f{R) (x R" 
with n = 1.1. Upper right plot: Same model as in the upper 
left plot but with C\ = C2 = po- Lower left plot: V{<j)) for 
w = 1/3 and 3m/2 = 1.1 in the case Ci = C2 = 0. This 
corresponds to f[R) (x ii" with n = 1.1. Lower right plot: 
Same model as in the lower left plot but with C\ = C2 = po- 
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FIG. 2: Left plot: V{(l)) for = and 3m/2 = 0.9 in the case 
Ci = C2 = 0. This corresponds to f{R) with the first term 
in (|lip proportional to ii" with n — 0.9. Right panel: Same 
model as in the left plot but with Ci — C2 ~ po- 



As long asm > 2/3(1 + w) (which corresponds to 
f{R) (X i?" with n > 1 in the case Ci = C2 = 0), the 
characteristic shape of the potential does no depend on 
the values of Ci , C2 , and w {w is in the range of phys- 
ical values < w < 1.). In any case, the scalar field 
starts at high absolute values and goes down its poten- 
tial to asymptotically freeze at = — 1. In the case 



m < 2/3(1 + w), the shape of the potential depends on 
the presence of non-zero Ci and C'2, as illustrated in Fig. 
[21 but the dynamics nevertheless drives (j) towards a con- 
stant value at late times. 



4 



DISCUSSION AND CONCLUSION 

We have shown here that exact power-law solutions in 
/(i?)-gravity can only exist for the very specific form of 
f{R) given in (fTT|) . If we ask for the theory to have the 
correct GR limit, then this f{R) simply reduces to i?". 
This makes R" gravity very special in the sense that it 
is the only f{R) model that allows for exact power-law 
solutions. Models with actions allowing for terms of the 



form R + f{R) for example, as studied in [15(, do not 
allow for exact power-law backgrounds unless we reduce 
to the GR background /(i?) = 0. 

We do not exclude the existence of cosmologically vi- 
able trajectories for a general f{R), but suggest that 
these trajectories may correspond to more complicated 
exact solutions that asymptotically scale like power-law 
solutions. 

We emphasize that one must make sure that such back- 
ground solutions exist before performing perturbation 
theory. It is not enough to simply perturb around an ex- 
act power-law background, since this exact background 
solution does not exist unless f{R) — i?". 

Furthermore, the above work suggests that in a dy- 
namical systems analysis of any f{R) theory other than 
i?" , one should not expect to find any equilibrium points 
corresponding to exact power-law solutions. 

We conclude that there is a qualitative difference be- 
tween i?" and any other /(i?)-model in the sense that 
i?" has exact power-law solutions even in the non- 
perturbative non-GR case (e.g. for n = 2), while any 
other f{R) can only allow these exact solutions in the 
GR-limit. Therefore, perturbations around these back- 
ground solutions should be carried out with caution. 
Moreover, we have shown that the singularity that ap- 
pears in other f{R) models is not present for the class of 
models derived from the requirement of a power-law ex- 
pansion, which is manifest in the scalar field analogy, in 
which the scalar field potential is a well behaved function 
of the value of the field cj). 
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